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ABSTRACT 

Seismic observations by the space-borne mission Kepler have shown that the core of red giant stars slows down while evolving, 
requiring an efficient physical mechanism to extract angular momentum from the inner layers. Current stellar evolution codes fail to 
reproduce the observed rotation rates by several orders of magnitude and instead predict a drastic spin-up of red giant cores. New 
efficient mechanisms of angular momentum transport are thus required. 

In this framework, our aim is to investigate the possibility that mixed modes extract angular momentum from the inner radiative 
regions of evolved low-mass stars. To this end, we consider the transformed Eulerian mean (TEM) formalism, which allows us to 
consider the combined effect of both the wave momentum flux in the mean angular momentum equation and the wave heat flux in 
the mean entropy equation as well as their interplay with the meridional circulation. In radiative layers of evolved low-mass stars, the 
quasi-adiabatic approximation, the limit of slow rotation, and the asymptotic regime can be applied for mixed modes and enable us 
to establish a prescription for the wave fluxes in the mean equations. The formalism is finally applied to a I.SMq benchmark model, 
representative of observed CoRoT and Kepler oscillating evolved stars. 

We show that the influence of the wave heat flux on the mean angular momentum is not negligible and that the overall effect of mixed 
modes is to extract angular momentum from the innermost region of the star. A quantitative and accurate estimate requires realistic 
values of mode amplitudes. This is provided in a companion paper. 
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1. Introduction 


Rotation has important consequences on stellar evolution. In 
particular, it induces meridional circulations, as well as shear 
and baroclinic instabilities, which contribute to the redistribu¬ 
tion of angular momentu m and to the mixing of chemical el¬ 
ement s (see fo r instance iTalonI I2008I: iMaederl l2009t iMathisI 
l2013t iPalaciosI l2013l for comprehensive reviews). Except for 
the Sun, the observational constraints remained sparse until 
the ad ve nt of the space-bor ne missions CoRoT ([ Baglin et"^ 


Bedding et al 


2006^^ Michel et al. 200 ^ and Kepler (iBorucki et al. l201(ii 


20 lot IChanlin et al.l l201 111 . These were mainly 


constraints on the efficien cy of the transport of chemicals (e.g., 
ICharbonnel & Zahnll2007h provided by observations of surface 
abundances. 

Based on seismic measu rements, iBeck et ^ (l2012ll and 
iDeheuvels et al.l (l2012l l2014t) brought stringent observational 
constraints on the rotation profiles in the innermost layers of sub¬ 
giant stars observed by Kepler. They concluded that the core of 
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subgiant stars spins up, while their envelope decelerates. Their 
core rotates about three to ten times faster than the envelope, 
depe nding on the evoluti onary state of the star. On the other 
hand, iMosser et al.l (l2012ll analysed a sample of about 300 red- 
giant stars observed by Kepler and found that, surprisingly, the 
mean core rotation rate decreases significantly during the red- 
giant phase. 


Current models of red-giant stars including angular momen¬ 
tum redistribution processes are unable to explain such low core 
rotation rates in subgiant and red giant stars. They are also 
unable to explain the deceleration of the core during the as ¬ 
cent of the red - giant branc h. Indeed, Eg genberger et al.l (l2012ll . 
iMaraues et al.l i2013h . and ICeillier et al.r(l2013 ) computed stel¬ 
lar models with transport of angular momentum by both merid¬ 
ional circulation and shear instabilities, concluding that these 
processes cannot explain the observed rotation profiles. The au¬ 
thors emphasised the need for a n additional phy s ical m echanism 
to do so. Along the same lines, ICantiello et al.l (l2014t) included 
the rotationally induced circulations and shear instabilities to¬ 
gether with the effect of a magnetic field generated through the 
Tayler-Spruit dynamo dSnrui ¥l999llT002h and reached the same 
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conclusion. However, a recent workbv rRiidiger et al.l(l2015h has 
shown that magneto-rotational instabilities of a toroidal mag¬ 
netic field could explain the angular mome ntum redistribution 
in sub giants and early red giants (see also iMaeder & Mevnej 
l2014l) . Interna l gravity wa vefl can also transport angular mo¬ 
mentum (e.g., iPressI [19811) and it has been demonstrated that 
they could explain the n early flat rotation profile in the inner ra¬ 
diative zone of the Sun (ICharbonnel & Tai^l2005h . While this 
process is st i ll to b e investigated for subgiant and red giant stars, 
[Fuller et^ (l2014l) have found that internal gravity waves are 
likely to couple the convective region and the upper radiative re¬ 
gion, but not the innermost layers. Hence, the question remains 
open, and we still need an efficient physical mechanism that ex¬ 
plains the slow-down in the core of evolved low-mass stars. 

Our aim in this paper is to investigate the influence of nor¬ 
mal modes (more precisely, mixed modes) on the core rota¬ 
tion of subgiant and red giant stars. Angular momentum can 
be transferred through energy exchanges between the oscilla¬ 
tions and the mean flow, which results in a modification of the 
rotation profile. Non-radial modes in subgiants and red giants 
have two interesting properties. Firstly, they have large ampli¬ 
tudes both in the outer layers and in the core of the star so that 
they are efficiently excited by turbulent convection while hav¬ 
ing a non-negligible effect on the innermost layers. Secondly, 
their amplitudes can be infer red, since they have been observed 
by CoRoT and Kepler (e.g.. Bedding et al.l 1201 ll: iMosser et alJ 
1201 It IChanlin & Miglioll2013l) . It is inossible. therefore, to pro¬ 
vide not only qualitative but also quantitative estimates of the 
angular momentum transport. 

The problem of the redistribution of angular momentum by 
progressive waves and normal modes has been considered for a 
long time fro m a th e oretic al point of view. For modes, the pio¬ 
neer work of lAndol (Il983h addressed the question of the inter¬ 
action of wave and rotation through wave momentum stresses 
in the mean angular momentum equation. Following this work, 
efforts have been made to theoretically address the problem of 
the redistribution of a ngular momentum in classical pulsators by 
unstable modes (e.g., Ando 198^ Lee_&_Sa^ 1993t lLeell2007l: 


iTownsend & MacDonaldl 2008t Townsendl 2014 ). Nevertheless, 
a quantitative estimate is difficult to make since it remains a 
challenge to predict the amplitude o f modes in classical pul 


sators (see, howev er, iLee et a l.ll2014t). In p arallel, [Press ([1981 ): 
Zalm et akl ([1997[) : [Kumar et al.l ( 199^ : [Mathis et al ( 20081 
20131) ha^ considered the interaction of the wave with the mean 


flow for progressive internal gravity waves (IGW) with a slightly 
different formalism. Nevertheless, all these works mainly fo¬ 
cused on the effect of waves on the mean angular momentum 
equation and neglected their effects on the mean energy equa¬ 
tion. 

In contrast, the interaction between waves/modes and merid¬ 
ional circulations has been extensively studied in geophysi¬ 
cal flows in the 60s and 70s in the context of middle at ¬ 
mosphere dynamics (e.g., [ Andrew s et al.[[l9^ [Holton[[l992[) . 
[Andrews & ]VlcIntvre[ ([19761 [1978bl) proposed the transformed 
Eulerian mean formalism (TEM) to account for the interaction 
between waves and the mean flow in a general way. That allowed 
them to demonstrate the non-acceleration theorem: if the waves 
are steady, conservative, and of small amplitudes, the mean flow 
is not accelerated. In the present series of papers, we adopt this 
formalism and adapt it to model angular momentum transport 


* Note that all along this paper, we will distinguish between progres¬ 
sive waves (hereafter mentioned by waves) and stationary waves (here¬ 
after mentioned by modes). 


by mixed modes in stars. In this paper, which is the first of the 
series, we establish the formalism and discuss numerical results 
for a benchmark model. In a companion paper (hereafter paper 
II) we consider realistic mode amplitudes to quantitatively esti¬ 
mate the efficiency of the angular momentum transport by modes 
along the evolution of low-mass stars. 

This paper is organised as follows: Section |2] introduces 
the mean flow equations. Section [3] discusses the transformed 
Eulerian mean (TEM) formalism and its particular form in the 
case of stellar shellular rotation. In Sect. |4] we model the wave 
fluxes for mixed modes. In Sect. |5] the rate of angular momen¬ 
tum transport is computed numerically for a benchmark model 
of an evolved 1.3Mo star, representative of CoRoT and Kepler 
observations. Section |6] is dedicated to discussions and conclu¬ 
sions. 


2. Preliminary remarks on the wave mean-flow 
equations 


In this section, we derive the equations that describe the inter¬ 
action between the mean flow and waves. We use an Eulerian 
averaging process to derive the equations of the mean flow. It 
allows us to discuss the coupling between waves and the merid¬ 
ional circulation as well as the role of the wave heat flux. 

The continuity, momentum, and energy equations in an iner¬ 
tial frame can be written as 

5+V-(pv) = 0 (1) 

Ot 

^^-(-V-(pvv) =-Vp-pVO + Z (2) 

Ot 

^ + v-(pv.) = e, (3) 

ot 

where p is the density, v the velocity field, p the pressure, s 
the specific entropy, O the gravitational potential, Z is a non¬ 
conservative mechanical forcing (e.g., turbulent dissipation), and 
Q represents heating or cooling terms. 

The azimuthal component of Eq. (|2]i can be written in a form 
that represents the conservation of specific angular momentum: 


diph) 

dt 


+ V ■ (p/rv) = - 


dp 50 
5 (^ ^ 50 


+ mX^, 


( 4 ) 


where the specific angular momentum is 5 = ur v^ = tir 
with Q representing the angular velocity, m = rsin 0, and v^ the 
azimuthal component of the velocity field. 

A given field can be decomposed into a mean part and a per¬ 
turbation. The perturbation is associated with linear non-radial 
waves. More precisely, for a given field A, we have 


A = A -(-A', 


( 5 ) 


where A is the Eulerian-mean azimuthal average 
- 1 

A = —Jo Ad0. (6) 

Introducing the decomposition given by Eq. (I5]l into Eqs. ([TJ, 
(|3]l, and (HJi and azimuthal averaging them, we obtain the follow¬ 
ing system: 

^ ■ (p vj = D (7) 

Ot 

dh 

p-^ +p (v± ■ Vi) 5 = -V_L ■ (czpv^Vi) + ra-Z^ H- TT (8) 

pf+p(r.'5’JT=-V,-(prtr) + e + A (9) 
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with 


"K = -p'y’±_ ■ Vj_/2 - 


S - -p'v^ ■ Vj_i - 




-_ \ dp's' 


dt 


[p's' Vi) 


(t3Tp'v;vi), 

( 10 ) 

( 11 ) 

( 12 ) 


and Vi, Vi ■ (), Vi are the gradient, divergence, and velocity 
vector in the meridional plane, respectively. They are defined by 


independent and, in the presence of wave forcing, the meridional 
circulation must ensure that h and s satisfy the baroclinic equa¬ 
tion. 

Consequently, both the wave momentum and wave heat 
fluxes must be considered to account for the effect of waves 
on the mean flow. To do so, we a dopt in the following section 
the TE M formalism introduced bv [Andrews & Mclntvrel (Il976[ 
Il978hh . 

3. Transformed Eulerian mean equations 

3.1. Formalism 


3 1 d 

Sr— + 

r dO 


dr 


and Vi = VrCr + Vgeg. 


(13) 


The terms involving the density perturbation (p') are gen¬ 
erally considered to be small (e.g., lAndol 119831 : lUnno et ^ 
Il989l) . This is particularly the case for low-frequency waves 
(cTji <sc N, where cr^ is the wave frequency and N the buoy- 
ancy frequency) where the a nelastic approximation applies (e.g., 
iDintrans & Rieutora 1200 Ih . Therefore, we neglect the terms 
p'v^, p's', and p'v'^ in Eqs. (fTOl i to (fT2l) . In addition, we use 
the Cowling approximation by neglecting the perturbation of the 
grayitational potential, so that Eqs. © to I© become 


^ + Vi-(pvi) = 0 (14) 

dt 

dh 

p-^ +P (Vi ■ Vi)fl = -Vi ■ (rrrpy^vl) (15) 

_ _ 

P^+p(Vi ■ Vi)s =-Vi ■ (ps'v^) H-2, (16) 


where the effect of wayes now appears through the waye mo¬ 
mentum flux in the mean angular momentum equation, Eq. (flSl l 
and the waye heat flux in the mean entropy equation, Eq. (fTfil l. 

Howeyer, the mean specific entropy s and mean specific an¬ 
gular momentum h are not independent. They are connected by 
the baroclinic equation (also known as the thermal wind balance 
equation), obtained by taking the curl of the hydrostatic equilib¬ 
rium equation: 


—2 w 

P Vi 


w 


-V0 


X VitjT + Vip X Vip = 0, 


(17) 


The waye heat flux R = i'v^ can be split into a component along 
an isentropic surfac e, the skew flu x, and a component perpendic¬ 
ular to it. Eollowing IVallisI (l2006h . 

R - (n X R) X n + (n ■ R) n , (18) 


with n - ViS/|ViS|. The diyergence of the skew flux can be 
rewritten as 


Vi ■ [{n xR)xn\ 



/ Vii X R 
I |Vi5p 
^ ViS X R 
IVi^P 



= V • ViS, 


(19) 


where we used the relation V ■ (a x Va) = Vo- • (V x a), yalid for 
any scalar a and yector a. 

Equation (fT^ shows that the skew flux behayes like an ad- 
yection by the yelocity v. The main motiyation underlying the 
TEM is thus to incorporate the adyectiye part of the waye heat 
flux into the mean meridional yelocity field (v). 

The residual meridional circulation is thus defined by 

pvi =pVi-H Vi x(p^e,^) , (20) 


where the stream function ij/ is defined, from Eq. (fT9l l. by 




Vii X R 
IViSp 


‘ ^(p 


1 

IViTp 



1 

r 



( 21 ) 


We note that the definition for t/r giyen i n Eq. (1211 1 is not unique 
(see lPlumb & Eerr^l2005l : I Valli J200^ for details). 

Inserting Eqs. (l20l i and (1211) into Eq. (fl4li to Eq. (fl6l l. we haye 


together with the equation of state 1 -s{p,'p). 

Therefore, we must solye the complete set of Eqs. (fill) to 
(flTl) to compute the effect of wayes on the mean flow, includ¬ 
ing both waye stresses in the mean angular momentum equa¬ 
tion, Eq. (fl5l) . and the waye heat flux in the mean entropy equa¬ 
tion, Eq. (UH). Howeyer, the effect of the waye heat flux is often 
neglected in the literature r elated to transport of angular momen- 
tum i n stars by vvayes (e.g., Press| [l981 [ Andojll^^ i Lee&S;a3 

a 


IZahn et al.ll 1 997l iKumar et al Jl 1 9991 iPantillon et al.ll2007l 


Mathis et al.ll2013l:lTownsendll2014). 

It is potentially misleading to oyerlook the effect of the waye 
heat flux. To illustrate this point, we consider the extreme case 
of a conseryatiye and stationary waye field. Under these assump¬ 
tions, the non-acceleration theorem states that the mean flow is 
not accelerated by wayes, although the waye momentum and 
heat fluxes defined in Eqs. (fl5l ) and (flfil ) do not yanish. In other 
words, the waye fluxes of heat and momentum produce a merid¬ 
ional circulation that cancels their tendency to affect the mean 
flow. Indeed, the waye momentum and waye heat fluxes are not 


^+V^-(pvi) = 0 (22) 

dh 

p-^+p {y[ ■'^±)Ti ^■ (pF) + TuX^ (23) 

n— 

P^+p(vl-V^)l=-V^-(pG) + e, (24) 

where the components of the yectors F and G are giyen by 


Fr 

Gr 


— iff dh 
v' + -_ 

* '■ r 39' 
il/ d's 


- — dh 

Fg - m Y.v„ - il/ — , 


dr 


(25) 

(26) 


Equation (flTl ) is left unmodified. 

Einally, when F and G are specified, the residual yelocity 
v)^ becomes part of the solution of Eqs. (l22l) to (l24l) together 
with Eq. (fTTI ). They are strictly equiyalent to Eqs. (fT4l ) - (fTfll ). but 
present seyeral conceptual and practical adyantages: 
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- The TEM equations enable distinguishing the advective and 
diffusive parts of the wave heat flux and incorporating the 
advective component i n the me an velocity field. 

- [Andrews & Mclnt\^ (Il978bh showed that V ■ (pF) and 
V ■ ipG) only depend on wave dissipation and non-steady 
terms. It is the non-acceleration theorem. Therefore, the 
TEM makes the adiabatic and non-adiabatic contributions of 
waves more explicit, the latter being able to modify the mean 
flow. 

- Einally, depending upon the considered problem, the TEM 
formalism allows us to gather in a single equation and in 
the single term F both the wave momentum and wave heat 
fluxes. It explicitly shows that wave momentum fluxes and 
wave heat fluxes do not influence the mean flow separately, 
but only in the combination given by F. Eor instance, in the 
quasi-geostrophic approximation G 0 so that the term F 
(the Eliassen-Palm flux) is the only wave contribution to the 
problem. As we will show in the companion paper, we are 
in the same situation for mixed modes in the limit of shel- 
lular rotation, in which the radial component of G can be 
neglected. 


3.2. Equation for the mean angular momentum: shellular 
rotation 


Eollowing the work o f lZah^ (Il992h: IV aeder & ZahnI (Il998h : 
iMathis & ZahnI (l2004l) (see also iMaedCT 20091 for an extensive 
discussion), shellular rotation has been an approximation widely 
used in ID stellar evolutionary codes. It is based on the as¬ 
sumption that t urbulence is highly anisotropic in stellar radia- 
tive zones (e.g.. lTalon & ZahtJl997tlMaedei)l2003HMathis et al.l 
I2OO4II : it is much stronger in the horizontal than in the vertical di¬ 
rection. Accordingly, efficient horizontal turbulent viscosity en¬ 
sures that the angular velocity is almost constant on isobars. We 
then consider the TEM equations in this framework. 

Given that the specific angular momentum is defined by h - 
uP' Q.{r, 6), we introduce the following decomposition 


The various physical quantities are expanded in Legendre 
polynomials 


UHr,6) = J] 

I 

Q (r, ff) - ^ ElcQ( (cos 6) 

t 

A(r,0) = 2AKr)EKcos0) , 


+ t dPf(cos 0) 

Uj (r)P ({cos 0)er + vj{r) -—— -eg 


d6 


(31) 

(32) 

(33) 


where Qe - Pe - 6(fi + (5f,2/5 (see lMathis & Zahnll2004h . and, 
for slow to moderate rotation, we assume that X c. (X). 

Eollowing these assumptions, the specific entropy is nearly 
constant on isobars and the stream function i/^ (Eq. 1211 1 becomes 
at dominant order 





V'i' 


(34) 


which is the same expression as for geostrophic flows (e.g., 
lHoltonlll992h . 

Introducing Eqs. OTI l - (l33l l into Eq. (l23l l. horizontal averag¬ 
ing, neglecting the time-variation of the density, and using the 
expression of the residual stream function given by Eq. (l34l i. we 
obtain, at dominant order. 


(P) 


d(r2Qo) 


1 d 


df dr 

with the fluxes defined by 

Tciic = --i (p) r-^lo u[ 


2 p. {r cue + ^hear + *ll'aves)j 


'Tshear = “ (p) Wr 


2 5O0 
~dP 


'll'aves ~ (p) ( ru 


—/d(5> 

c' I - 


dr 


(35) 


(36) 

(37) 

(38) 


Q(r, ff) = Qo (r) -1- D. (r, 0) 


(27) 


with 


Q() - 


sin^ 0 Q(r, 0) d0 
sin^ 0 d0 


3 r . 

4 Jo 


sim 0 Q(r, 0) d0, (28) 


where, with the shellular approximation, Qq » G. The scalar 
quantities are developed as 


The turbulent stresses are described by an anisotropic eddy vis¬ 
cosity denoted Vv and vi, for the vertical and horizontal compo- 
nent, respectively (see IZahiilll 9921 iMathis & Zahnll2004 for de¬ 
tails). We also note that we introduced the Lagrangian derivative 
d/dt^d/dt-hrd/dr. 

Under the assumptions used for Eq. (l35l l. we obtain for the 
entropy equation 


with 


(P) 


d<5) 

df 


11 . 

r“ dr 


(r^S}-i-(Q) 


(39) 


X(r,0)^(X)(r)+X(r,0) , (29) 


.s = (p> s' v;, 


(40) 


with 


(A> = 


sin 0 X(r, 0) d0 
sin 0 d0 


_ 1 r 

“ 2 Jo 


sin0X(r,0)d0, 


(30) 


and the velocity held equals to (r, 0) = rCr -i- (r, 0), where 

f comes from the contraction or dilatation of the star and is 
the residual meridional circulation. 


and 


{T}{Q}^{ps}+^_j^l^r^(x)^Y (41) 

where T is the temperature, e is the nuclear energy genera- 
tio n rate, and y is the t hermal co nductivity (see, e.g.. S ect. 6 
of IMathis & Z^ 12004 and also iMaeder & Zahnlll998h . This 
expresses energy conservation on a level surface. To obtain 
Eq. ( l39l l. we used Eq. (l34l i and assumed that the horizontal 
derivative of entropy is negligible at dominant order. 
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3.3. Discussion 


4.1. Simplifying approximations 


As mentioned in Sect. |2] transport of angular momentum by 


ing the effect of t he wave heat flux (e.g. , 
iLee &Saiolfl^ IZahn et al.l[T^ 
Townsend Fo llowing the work of 


PressI 1981 

-- 

Andol 

Kumar et al. 

|l99^ 

Bretherton ( 

19691). 


wave 


contribution in the mo mentum equation, equal to 2 cos 6 Qq v'fg 
in our notation (see also lLeell200^lMathisll2009t) . This term cor¬ 
responds to a Lagrangian correction that can be understood as 
follows: on the typical time-scale of the momentum transport by 
waves (the modal period, for example), rotation advects a wave- 
related momentum flux. This term arises naturally in Lagrangian 
mean wave-flow theories su ch as the generalized Lagran gian 
mea n (GLM) formulation by Andrews & Mclnt\^ (Il978al) (see 


also lGrimshawll 1984t iLeelbOl 3l) . 

For linear waves, one can show that this Lagrangian term 
can be recovered from the TEM equations if one assumes that 
the Lagrangian perturbation of entropy vanishes {i.e. ds - 0). 
The wave flux in the mean angular momentum equation (Eq.ILSTi 
is 


v^v;-H2cos6>Qov^s'^^^j . (42) 

We can then express the entropy Eulerian perturbation in 
terms of the Lagrangian perturbation at first order: bi = -H 
frd('S) /dr, where is the wave displacement. Using Eg. (l45l l 
and the expression for the wave velocity v' = i {ctr + mQo) ^ - 
m (^ • VQq) (where ctr is the modal frequency), we obtain 


We restrict ourselves to the case of strictly shellular rotation (see 
Sect. I3.21 i. Moreover, we focus on the inner radiative regions 
since our aim in this series of papers is to investigate the trans¬ 
port of angular momentum in the core of evolved low-mass stars. 
These restrictions allow us to use several approximations to de¬ 
scribe the wave field: 


1 . 


2 . 


3. 


The quasi-adiabatic approach: It consists of neglecting the 
difference between adiabatic and non-adiabatic eigenfunc¬ 
tions in the full wave equations. This approximation is valid 
when the local thermal time-scale is much longer than the 
modal period. This is the case in the radiative region of 
evolved low-mass stars. 

The low-rotation limit: We assume that the modal period is 
much shorter than the rotation period. This is justified by 
recent inferences of the rota tion in the core of subgiants 
( Deheuvels et al.l 120121 l2014l) and red giants (iMosser et al.l 
12012 ) using seismic constraints from Kepler. 

The asymptotic limit: Finally, we use an asymptotic de¬ 
scription for gravi ty modes (e.g.. lDziembowski et al.ll200ll : 
iGodart et al.ll2009li . which is valid for mixed modes in the 
inner radiative re gion of subgiants and red giants (e.g., 
iGounil et al.ll2013h . 


4.2. Derivation of the wave fluxes 

We assume that the perturbations consist of non-radial oscilla¬ 
tions such that, for any perturbed variable /', we have 


(43) 

Therefore, if the entropy Lagrangian perturbation vanishes (i.e. 
6s = 0), the term (l42l i becomes 


v^v;H-2cos6»Qov'^e, (44) 

which is the expression given in iPantillon et al.l (l2007l) as 
claimed. Fina lly. Eg. (l44b shows tha t the Lagrangian correction 
introduced by IPantillon et al.l (l2007l) ; iLeel (l2008h : lMathisl (l2009l) 
is recovered from the TEM equations. In this framework, it can 
be understood as originating from the wave heat source term in 
the mean entropy equation. 

However, the contribution related to the entropy Lagrangian 
perturbation is not negligible (see the results presented in 
Sect. l5.2l below). Therefore, it is not enough to consider the wave 
flux in the mean angular momentum given by Eq. (l44l i. 

4. Modelling the wave fluxes: the case of mixed 
modes in evolved low-mass stars 

In this section, we focus on the transport of angular momentum 
by mixed normal modes in the radiative region of evolved low- 
mass stars (i.e. subgiants and red giants). These modes have a 
dual nature; they behave as acoustic modes in the upper layers 
and as gravity modes in the inner layers. Therefore, they are de¬ 
tectable at the stellar surface while providing information on the 
innermost regions. Th ey have been detected in several thousands 
of evolved stars (e.g.. llBedding et aT1l201 ll : [Mosser et al.ll2012l 
l2014l) and have provided a wealth o f information on the core of 
red giants (e.g.. iMosser et al.]l2012l) . These modes must not be 
confused with progressive gravity waves, often called internal 
gravity waves in the literature. 


f =aK {/'(r, 6) , (45) 


where Tie denotes the real part, a is the amplitude, and ctr is the 
real frequency (so that non-adiabatic effects are neglected). The 
amplitude is considered as statistically constant in time since we 
deal with solar-like oscillations resulting from a balance between 
mode driving and damping (see lSamadill201 ll for details). 

We now consider the equations governing non-radial and 
non-adiabatic oscill ations in the prese nce of rotation (e.g., Eqs. 
36.12 to 36.15 of lUnno et al.l Il98^ . Eor sake of simplicity, 
tildes, bars, and brackets will be dropped from now on unless 
necessary. The equations are 


Dvl 2/r 1„ , p' 

Df UT m p 

Dip' p' > A 

Df m ^ 


where Fi = {d\np/d\np)^, pj = - {d\np/d\nT) and 


(46) 

(47) 

(48) 

(49) 


D 

Df 




A = Vlnp-Vln p , 

Ti 


(50) 


so that = p 'Vp ■ A. 
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Given the assumption that the mode amplitude is statistically 
constant in time, Eqs. (l46l l to (l49l) reduce to 

= --V_l/7' + ^V_lP (51) 

^ p 

1 im d' 

i{(TR + mQo) v; + -vl ■ V^/r =-^ (52) 

^ VJ TU p 

/(ctr + mQo) H-pr—) =-vl • A (53) 

\P Tip Cp) 

iio-R + mno)p' + V_L ■ (p\'J + —vj, = 0. (54) 

ZJ 

Using Eqs. (EB to (l54l) and the assumptions detailed in 
Sect. 14. II and after stirring vigorously as detailed in AppendixiAl 
we obtain for a mode of a given angular degree { and azimuthal 
order m an expression of the form 



j__5 

f.2 Qp 



(55) 


a 


2 

{,m 



5“ (r^Qo) 
dr^ 


+ &; 


5(r^Qo) 


+ Cl r^Qo -H md-Dj 


Fig. 1. Rotation rate normalised by its central value (solid line) 
and buoyancy frequency normalised by its maximum value 
(dashed dots lines) versus normalised radius for the model de¬ 
scribed in Sect. 15 .1 1 (logarithmic scale in both axes). 


where & - ctr + mQo and the coefficients JV'i ,S”i ,Ci ,Tfi are 
given by Eqs. (IA.25I I to (IA.28I) . respectively. The amplitude 
corresponds to the amplitude of a mode of angular degree (, and 
azimuthal order m. 

The derivation of the wave fluxes in Eq. ( l39l l is more direct. 
Using the asymptotic expression for the Lagrangian perturbation 
of entropy given by Eq. (IA.22I) . we obtain 


1 d , 

C dr V 


d 


di, 




dr 


(56) 


where 


a - - 


Anr^pT \ V 

{{{+1) (N^ 




dr 


\^R 


(57) 

(58) 


with L being the luminosity, N the buoyancy frequency, V and 
Vad the actual and adiabatic temperature gradients, respectively. 


5. Case of an evolved 1.3 M© star 

The aim of this section is to qualitatively discuss the effect of 
mixed modes on the mean angular momentum profile within the 
TEM framework. To this end, we consider a stellar model for 
an evolved I.SMq star that is representative of the oscillating 
red giant stars observed by CoRoT and Kepler. A quantitative 
estimate as well as a discussion of the effect of mixed modes 
throughout the evolution is given in paper II. 

5. 1. Computing the benchmark model 

The equilibriu m model was compute d with the stellar evolution 
code CESTAM dMaraues et aklTOl^ . The atmosphere was com¬ 
puted assuming a grey E ddington approxima tion. Convection 
was included according to ICanuto et al.l (Il996ll . with a mixing- 
length parameter a = 0.67. The initial chemical composition 
was the solar composition of lAsplund et al.l (l2005l) . with a he¬ 
lium mass fraction Y = 0.261 and metallicity Z = 0.0138, and 


Table 1. Characteristics of the models computed with rotation 
and without rotation for comparison. 


Model 

R/Re 

LILq 

Teff (K) 

Age (Myr) 

Rotating 

3.558 

6.830 

4952 

4716 

Non-rotating 

3.555 

6.784 

4945 

4590 


diffusion was not inc luded. We used t he OPAL equation of stat e 
(iRogers et al.l 1 19961) and opacities (llglesias & Roger j 199^ 
compl emented, at T < 10"^ K, by the Alexander & EergusonI 
( 19941) opacities. We adopted the NACRE nuclear reaction rates 
( Angulo et al.lll99^ exc ept for the + p re action, where we 
used the reaction rates of llmbriani et al.l (l2004l) . The model was 
computed without rotation. 

To obtain a rotation profile as realistic as possible, we com¬ 
puted an eq uivalent model with rotation (using the procedure 
described in iMaraues et al.ll2013l) and rescaled the rotation pro¬ 
file to gain a core rotation rate Qc/(27r) = 1 pHz. This value is 
consistent with the mean central rotation rate of evolved low- 
mass stars derived from Kep ler observations (see iMosser et al.l 
120121 iDeheuvels et al.ll20l4) . The resulting rotation profile and 
buoyancy frequency are shown in Eig. [1] Consistently with the 
slow rotation limit we consider, this procedure is equivalent to 
assuming that rotation has a negligible effect on the star equi¬ 
librium structure. The equivalent model including rotation was 
computed using the same input physics and parameters as the 
model without rotation. It has a radius as close as possible to the 
radius of the model without rotation and is in the same evolu¬ 
tionary phase. The main characteristics of the models are given 
in Table [1] 

Einally, we used the ADIPLS co de for adiabatic oscillations 
(IChristensen-Dalsgaardll20()^ 1201 ll) to compute eigenfunctions 
and eigenfrequencies. The effect of rotation was also neglected 
in computing the eigenfunctions and eigenfrequencies. 


5.2. Effect of prograde and retrograde mixed modes 

We will show in Paper II that the wave fluxes in the mean energy 
equation, Eq. (l39l) . is negligible. Therefore, within the TEM for- 
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malism, the contribution of the mixed modes to the mean flow 
consists entirely of the wave flux in the mean angular momen¬ 
tum equation (see Eg. [15]). Thus, for ease of notation, we define 
the rate of temporal variation of the mean specific angular mo¬ 
mentum induced by mixed modes by 

j ^ , (59) 


where we recall that iFwaves is given by Eq. (l38l l. Note that J is 
normalised so that for a given mode of angular degree { and 
azimuthal order m one has IR\ = 1 at the photosphere. 

Computing j as described in Sect. 14.21 leads to the conclu¬ 
sion that prograde modes {m < 0) extract angular momentum 
in the regions near the maximum of the buoyancy frequency 
and the maximum of the rotation rate, therefore slowing down 
the core. Conversely, the retrograde modes (m > 0) tend to 
spin up the core. This is illustrated in Eig. |2] (top and mid¬ 
dle panels) for the sectoral modes ( - 2, m - {-2, -i-2) at 
the frequency v« = cri(l{2n) - 300.9//Hz, but this conclusion 
is valid for the other modes as well. When we consider un¬ 
stable modes, the opposite situation is found, th at is, the ret¬ 
rograde modes ex t ract angular momentum ( e.g., lAnd^ 198^ 
Lee & Sai^ j_993^ Townsend & MacDonald! l2008r Townsend 


2014HLee et al.l2014l) . In this case, the wave stresses in the mean 


angular momentum equation are dominated by transient terms. 

We have also found that, despite an important rotation gra¬ 
dient in the hydrogen burning shell, the dominant contribution 
to the angular momentum transport comes from the term 2)^" so 
that, from Eq. (i55l l and for a mode of a given angular and az¬ 
imuthal degree, we can write 


j ^ m& , (60) 

where 29'^' is given by Eq. (IA.28b . 

To proceed, we considered the sum of the contributions of 
a prograde and a retrograde mixed modes to the mean angular 
momentum profile. In the absence of rotation, both prograde and 
retrograde modes transport angular momentum, but the sum van¬ 
ishes so that there is no net transport. Including rotation, the net 
transport of angular momentum is a small residual (see Eig. 
bottom panel). The asymmetry between prograde and retrograde 
modes arises from the Doppler shift terms in the mode frequen¬ 
cies {i.e., & - cTj^ + More precisely, the dominant term of 

Eq. (i60l l comes from the first term ofin Eq. (IA.28b . This term 
originates in the second term of Eq. (i38b . thus from the wave heat 
flux. To illustrate it, we added a prograde and retrograde mode, 
expand in flo/o"/?, and only kept the leading term to obtain 


j{i, -\m\) + J(£, \m\) ^ 2\m\Wral\m\\^r''"f 




( 61 ) 


which is negative, since the term a (see Eq. [57l i is nega¬ 
tive. The sum of the contributions of prograde and retrograde 
mixed modes thus implies an extraction of angular momentum. 
Nonetheless, in the regions where the contribution of the pro¬ 
grade mode has nodes, the contribution of the retrograde mode 
dominates due to the second and third terms of 2?^' in Eq. (IA.28b . 
It therefore results in a complex situation where there is a spa¬ 
tial alternation between regions where mixed modes decrease 
and increase the mean momentum, as shown in Eig. |2] (bottom 
panel). However, before reaching a conclusion, we have to con¬ 
sider the collective effect of mixed modes. 



r/R 




r/R 


Fig. 2. Rate of angular momentum transport by mixed modes {j 
as defined by Eq. l59l l computed as described in Sect. 14.21 for 
the \3Mq model described in Sect. O versus the radius nor¬ 
malised by the radius of the radiative region, j is normalised 
so that Iff’™//?I = 1 at the photosphere. The colour code is as 
follows; the solid red line corresponds to 7 < 0 (thus to a spin 
down) and the solid blue line to 7 > 0 (thus to a spin up). The top 
panel exhibits the result for the prograde mode € - -m = 2 of 
frequency v« = 300,9 //Hz, the middle panel corresponds to the 
retrograde mode i - m -2, and the bottom panel to summation 
of 7 for the prograde and retrograde modes. We only show the 
inner radiative region, for which our formalism is valid. 
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r/R 


Fig. 3. Top panel: Mode inertia versus mode frequencies for an¬ 
gular degrees £ - 0 and £ - 2. Mode inertias are computed by 
integrating the squared eigendisplacement over the stellar mass. 
Bottom panel: Same as Fig. ^ except that j is computed by 
considering the total contribution of the £ - 2 modes for the 
selected range displayed in red in the top panel. 


5.3. Collective influence of mixed mode 

Oscillation spectra of low-mass evolved stars exhibit quasi- 
periodic patterns with periods that correspond to the large sep¬ 
aration as depicted in Fig. [3] (top panel). Each pattern is made 
of p-dominated mixed modes separated by several g-dominated 
modes. It is therefore instructive to consider the collective effect 
of mixed modes on the mean angular momentum. 

We selected a range of modes, shown in red in Fig. [3] (top 
panel), and computed the total contribution to J by summing the 
contribution of all modes belonging to this range, including the 
different values of the azimuthal order m for each mode. Given 
the spatial shift between the consecutive mixed modes, the sum¬ 
mation induces a smoothing of J, as shown by Fig. [3] (bottom 
panel). The collective effect of mixed modes is to decrease the 
mean angular momentum and thus to slow down the core rota¬ 
tion of the star. 

Nevertheless, we emphasise that to realistically compute the 
effect of the modes on the mean angular momentum profile, we 
need to consider all modes with realistic individual amplitudes 
This is the scope of paper II. 


6. Conclusion 

Our aim here was to establish a formalism allowing us to calcu¬ 
late the transport of angular momentum by mixed modes in the 
inner radiative region of evolved low-mass stars. We thus con¬ 
sidered the transformed Eulerian mean (TEM) formalism that 
allowed us to take into account the combined effect of both the 
wave momentum flux in the mean angular momentum equation 
and the wave heat flux in the mean energy equation. Indeed, the 
effect of the wave heat flux on the mean flow is generally ne¬ 
glected in the context of angular momentum transport by waves 
in stars. Here, the TEM allows us to include its impact on the 
mean angular momentum in a relatively simple way. 

Given the important contraction of the core and expansion 
of the envelope of low-mass stars after the main sequence, a 
strong radial rotation gradient appears. This was taken into ac¬ 
count by adopting the shelullar approximation. While a two- 
dime nsional treatment is certainly des irable and more accurate 
(e.g.. lEsr)inosa Lara & Rieutordll201^ . this approximation pro¬ 
vides a first approach of the problem while keeping it tractable. 

Subsequently, the wave field has been modelled so as to ac¬ 
count for the transport of angular momentum by mixed modes. 
We limited our investigation to the inner radiative regions. Thus, 
several simplifications are justified; the low-frequency limit (i.e. 
cr“ «: {N~,S^}), the quasi-adiabatic approach (the local ther¬ 
mal time-scale is higher than the modal period), and the low- 
rotation limit as suggested by recent astero seismic observations 
(iDeheuvels et al.ll2012HMosser et al.ll2012h . It allowed us to ob¬ 
tain an explicit expression of the wave-related terms appearing 
in the mean angular momentum equation. 

Then, the formalism was applied to a benchmark model of a 
1 .3Mq evolved star. We found that the influence of the wave heat 
flux on the mean angular momentum is not negligible and even 
dominant when considering the sum of prograde and retrograde 
modes. In addition, we showed that the overall effect of mixed 
modes is to extract angular momentum in the innermost radiative 
layers. 

Nevertheless, more realistic estimates of mode amplitudes 
are required for a more quantitative determination of the amount 
of angular momentum transported by mixed modes. The second 
paper of this series is dedicated to this. 
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Appendix A: Derivation of the wave flux in the 
mean anguiar momentum equation (Eq.|35) 

We assume that the perturbations consist in non-radial oscilla¬ 
tions and that they behave linearly as exp [/(erf -h ot(^)], where 
cr is the modal frequency, t the time, and m the azimuthal de¬ 
gree. Consequently, for a given quantity /, its perturbation can 
be written 

f [fir, 0) , (A. 1) 

where a is the amplitude, the real part of the frequency cr, 
and 'Re stands for the real value. From Eq. (lA.lb one immediately 
deduces a relation for squared quantities that reads 

iJr -a^'Re[fr] ■ (A.2) 

One can then express the Eulerian perturbation of the entropy 
in terms of the Lagrangian one, using 6s - s'+^r d {s) /dr, where 
is the radial component of the wave displacement. It allows 
us to rewrite the divergence of the wave flux in the momentum 
equation (see Eq. [39l l such as 


.2 ('■^'^waves) - , [{fr^TU (fi) 5l)} (A.3) 


1 d 


'■ dr 


2r^ dr 


with 


= v/v? - 2cos0Qo 


-A -V* — 


-1 


(A.4) 


In the following, our objective is to provide a general expres¬ 
sion of Eq. (IA.3b using the full set of equations governing waves. 
Note that in the following, the tilde will be omitted for the per¬ 
turbation as well as overbar and brackets for the unperturbed 
quantities. 


A.1. Using the full wave equations for expressing Eq. lA.4tl 

In this subsection, we aim to express the correlation products 
appearing in Eq. (IA.4b using the set of Eq. (fSTT) to Eq. (l54l i. 


A.1.1. Derivation of the first term of Eq. IA.4b 

Eirst, Eqs. (l52l i and (l53]) are inserted into (l54l) . and an integration 
by part is performed. Second, Eqs. dSB to (l5^ are used as well 
as the relation 


■WF'' ''' ' 

. (A.5) 
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Finally, taking the real part and considering that the imaginary 
part of the frequency vanishes (ct/ = 0) since we consider steady 
waves (see Sect. lO l. one gets 


2ut \ m^Y 1 p &■ 



mp 

[(vlp-'V.p)vl--/l] 

- 

2 


1 ) 


Ai1. ■ 


[(TV 


well as Eqs. (ISTT) to to get 

p' iV/i /V ■ (pv^) 6s 

■p— =-^ ;-r +P7- — 

1 1/7 CT per Tucr Cp 


and 


pa-- 


sin^ 6 dp^ A dp 


cl dr dr 


dr 


P 

(vl ■ A + vl ■ Vi) (vl* ■ Vi/r) 

pm 

ds 


U r. • wdpl 

2 

& 

2 

^(p ^Pt 

+ 

/ 

2pQo smd-T — 

mc^ or 


(A. 10) 


(A.ll) 


, 5prsin0cos0Qo , dp 6s 

- -^l - '^^-P^TrTp- 


mp 


(vl -p ^'^±p)pTds*ICp 


(A.6) 


Finally, inserting Eqs. (lA.llI) and (lA.lOb into Eq. ( IA.9b . in¬ 
tegrating over the solid angle and taking the real part gives 


where we introduced the notation & - cr + inYlo, pr - 
-(d\vipld\TiT)p, and the subscripts I and R denote the imag¬ 
inary and real part, resp ectively. Note that this expression corre- 
sp onds to Eq. (3 6.16) of lUnno et al.ld 19891) and was first derived 
bv lAnd^dlgO . 

To go further, we assume shellular rotation and that the equi¬ 
librium state is not deformed by rotation. Moreover, integration 
over the solid angle is performed so that it finally gives 


15/, , pr 5(r2Qo) 




2r'^ dr 
pr^Yli 


-{sm^ecos6\'gv'p 


* 2mcl dr 

mpQ.0 


-CT 


(sin^ev^v;)^ 


+ 

^ cr 


mp ^ p(Ar- l)flo 

■(sm5v,v,)^-- 


2rd- dr 

prYlo I , dv 

—;— (cos0sm0v^ 

prpT 


p 5(r2Qo) / . , dv'; 
dr I, ■ ^ dr 


, . . I ■ A , mprdp I ,6s' 

+ —^cr(sin0—Vrf, ) (v 


Cp Ij 26- dr \ Cp Ip 


A.1.2. Derivation of the second term of Eq. fA.4b 
We start from Eq. d52l i. which immediately gives 


2Q() cos 0 Vg = 


imp 

mp 


■ icrv. 


/sind) 


(r^Qo) 

I ^ J 

' dr 


(2f2or^nrp cos 0Vg^'*) = -i, 
dv 


dp' 

r^ dr '' dr 




d^ (r^Ylo) 


- p sim 0 


d{r^Ylo) 


dr 




r'^ dr 


!sin6»p(Tv^f; 


51np 2crR m 


5 In r 6- r& dr 


1 d 


2ri (2Qor-"uTpcos Ov'^C)p = 
d 


Ar -H 3 -b 


r& dr 

mpcl I 1 d(sm0vg) 
2r^& \sin0 dd 


!■ a A '"Fr dp j 6s , 
^sin6»V0V^)^ - -TT—— V, 


2d- dr \cp 
6s 


m 

+ TTPrEip 


d{r^C) 


-P dr I j 


Pf ^ . 

-H —cr(sm 0 


dr dv'. 
dr dr 


mWj 1 ,d{rH'i) 


j 2r^& \sin^ dr jj 


(cos 0 v'gV';)^ 

(cosdsind v'vg*)^ 


pQo 1 5 In p 
& Yi dinr 


(cos d sin 0 VgV*) 


/ ■ 


(A. 12) 


A.1.3. Derivation of the third term of Eq. dA.4b 
The derivation is similar to the previous subsection, it gives 
Id 


/ 

(A.7) 


d (r^Qo) 


1 _ ^ dp 
pa-^ dr 


Wr^R 


P. 

2^ 


Ar + 3 + 


r& dr 


2 


1 d (r^Qo) d^ (r^Qo) 

(Ar + 2)-+ ^- L 

r or 


(sin^dv;^:)^ 


dr2 


v' (A.8) +p 


pdlr^Ylo) I I dr dV 


Equation Eq. dA.8b is further multiplied by r^pm^; and radial 
derivation is performed to give 


pr . 


, ,.AA\\ . mAYipl ,d(r2^:) 

\'''' dr 




2r'^& 


Vi ■ (pVi) 


d(r2f:) 


2 2 
'« pc. 


dr Ip 2r^& \sind 


w- 


d(r2^:) 


dr 


rfin dp , , 

— ^ (sindcosdvgf,)^ , 


where we have introduced the notations 


6 = 1^1 ds and v^ = fd-f.. 


(A. 13) 


(A. 14) 


(A.9) 


Now one needs to express both the perturbation of pressure and 
its radial derivative. To this end, we use Eq. (l53T l and Eq. (l54l i as 


A.2. The case of low rotation and frequency limit 

As described in the main text (see Sect. 14.2b we restrict ourselves 
in the limit of low rotation (i.e. ctr » Qq) and low frequency 
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(i.e. ctr «c N, where N is the buoyancy frequency). Those ap¬ 
proximations will permit us to derive tractable expressions for 
Eqs. (IA.7b . (IA.12b . and (IA.13b . To this end, the first step consists 
in considering the wave equations without rotation and introduc¬ 
ing them into Eqs. (IA.7b . (IA.12b . and (IA.13b . This is equivalent 
in considering a first-order development in term of rotation for 
Ea.dAAl). 

Projection onto the spherical harmonics is thus performed in 
the limit of low rotation. Eor one normal mode of a given m and 
{, the eigen-displacement and velocity can be decomposed such 
as 


? - ^{,'n + S/j ^{,m , 

where (lRieutordll987h 


^i,m ~ Tf ^ r 5 and 


dT” 


eg + - 


1 dV! 


do sin 0 d(p 


so that 


yni _ ^c,m _c_ w _ 

Sr - Sr ^ Se - Sa ’ S0 “ S/, 


do 


and 


dY"; 


1 dY'l' 
sin 6* d(p 

1 dY”' 


v; = !«'; Yf iw’” iV” -— 

" ^ f ’ e f do ^ sine d4> 


(A. 15) 

(A. 16) 
(A. 17) 

(A.18) 

(A. 19) 


where we used v"' = 

In the non-rotating limit and providing the decomposition 
given by Eq. (IA.15b . using Eqs. (fSTI) to (l54l i. v™ can be expressed 
as a function of di™ and only through the relation 


V 


m 

e 


= r@ 


O'rPt - 

On 


with 




j-2 ^-2 

' Or 



1 d(r2M'; 

^ u'f d In p 

(A. 20 ) 

y.2 dr 

1 

0 - 


2 


,and S ^ 

= 7(7+1)^. 
r^ 

(A. 21 ) 


To go further, in the asymptotic, quasi-adiabatic limit (ctr <sc 
N), 6s”' is a function of (o r equivalently so that v™ de¬ 
pends only on m™. It reads (see lGodart et al.ll2009l for details) 


and 




(A.22) 


dr r '' 


(A.23) 


where L is the luminosity, T the temperature, p the density, V 
the temperature gradient, and Vad the adiabatic temperature gra¬ 
dient. 

Einally, after projection onto the spherical harmonics and 
using Eq. (IA.20b . Eq. (IA.22b . and Eq. (IA.23b . one obtains for a 
mode of a given angular degree {€) and azimuthal degree (m) an 
expression of the form 


-i^rV 

Or, ' y w: 

or 


(A.24) 


fjiif 

JXp 


5“(r^Qo) 5(r^Qo) 


dr^ 


+ B’ 


dr 


+ C'i' r"Qo ■ 


• md" lyi 


with 




(A.25) 


47rSr,>f 


■C,mn 


. l^(A + 2^, + Aln[ry 7^, 


(A.26) 


r^a-Rcr 


(Tr / , K()\ 


2c? 


47rCf,,n _ 

pkj\^/r 

+ 

cr 


O-R 1 dp 

- — ■p — ■^@K2 (fo - Q'6) 
cr V\p ar 

(Ar 1) d 


(A.27) 


dr 


and 


47r £)/•.„, 
pkj\^;y 

- 


p&^ dr 


(A.28) 




1 dlnp 
n dr 


+ 


A -H - -H 2^1 -H ^ In {rQo-R^z) + :r ^^ 2 ) 
r dr dr 


where for Eq. (IA.28b we use the approximation Q ' = -i{{ H- 1) 
obtained in the limit 5 ^ » cr^, and 

,/v._iW^r' 


Anr'^pT \ V 


^ Pt Ids 

fo = a— — 
Cn \dr 


■t,r. 


Cl - 


akj 2 

^2 = — ; 6 = - + A + p . 

(Tr r 1 1 dr 


dr 

1 dlnp 


(A.29) 


as well as 


+ ({ +If {J7f 

K2^e{Jtif-(e+i){J7f- 


(A.30) 


Einally, we note that the normalization condition for spheri¬ 
cal harmonics is (Yf'Y’f'j — 6m,m’, where 6 is the Kronecker 

symbol. Note also that the following relation have been used 


nrvc v”’ — r'" V'" - 1 - r'«y»' 
Losor^ - ■'r+i-'r+i■'r''r-1, 

dY’” 


(A.31) 

(A.32) 


d^Y'y cos 
-I- 


OdY'I 


-H 


1 d'^Y’f 


dO^ sin 0 do sin^ 0 d(p^ 
where 

J f — 


^-{{(+l)Y”\ (A.33) 


(4(’2 - 1) 


1/2 


(A.34) 


if ^ > \m\, and 7? = 0 otherwise. 


































































